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Abstract 

Using the relation between the structural and the electronic properties of hon- 
eycomb, we study the hidden SU (3) symmetry of the graphene monolayer and 
exhibit the link with its electronic properties. We show that the conservation law 
of incoming and outgoing electronic momenta at each site of graphene is solved 
in terms of SU(3) representations; and the Fourier waves 4> (k x , k y ) of the hopping 
electron may be classified by SU(3) highest weight multiplets </> Pig (£)- ^ is also 
shown that the phases arctan y- of the waves are quantized as with p, q 

positive integers. Other features are also discussed. 

1 Introduction 

Graphene is a system of carbon atoms in the sp 2 hybridization forming a 2D honeycomb 
lattice. This is a planar system made of two triangular sublattices A and B and constitute 
the building block of graphite. Since its experimental realization in 2004, the study of 
the electronic properties of graphene with and without external fields has been a big 
subject of interest some of its main physical aspects were recently reviewed in [5] 

and refs therein. This big interest into graphene and its derivatives is because they offer 
a real alternative for silicon based technology and bring together issues from condensed 
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matter and high energy physics [6]- [12] allowing a better understanding of the electronic 
band structure as well as their special properties. 

In this paper, we focus on an unexplored issue of 2D graphene by studying the link 
between specific electronic properties and a class of hidden symmetries of the honeycomb. 
These symmetries allow to get more insight in the transport property of the electronic 
wave modes and may be used to approach the defects and boundaries [13] . The existence 
of these hidden symmetries; in particular of the remarkable hidden SU (3) invariance 
considered in this study, may be motivated from several views. For instance from the 
structure of the first nearest carbon neighbors like for the typical (A -Bi), (A -B 2 ) , 
(A -B 3 ) see fig(EE]) for illustration and more details. These are basic patterns generating 




Figure 1: Sublattices A and B of the honeycomb. A-type carbons are given by red balls and 
B-type atoms by blue ones. Each carbon has 3 first nearest neighbors and 6 second nearest 
ones. 

three SU (2) symmetries contained in the hidden SU(3) invariance of honeycomb. The 
(A-B) patterns transform in the isospin | representations and describe the electronic 
wave doublets ± i = [a (r) , b (r)] interpreted as quasi-relativistic 2D spinors in the 
nearby of the Dirac points [5]. The SU (3) hidden symmetry of honeycomb is also 
encoded in the second nearest neighbors ((Ao-Ai)) and ((B -Bi)), i = 1,...,6, which 
capture information on its adjoint representation where the six ((Ao-Ai)) (and similarly 
for ((B -Bi))) are associated with the six roots of SU (3). In addition to above mentioned 
properties, hidden symmetries of graphene are also present in the framework of the tight 
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binding model with hamiltonian, 

H= -h J2 i2a(Ti)tf(Ti+6 n ) + hc 

n€A n=l 

-h £ «(r i )at(r j ) + 6(r i )6t( rj ) , 

where the fermionic creation and annihilation operators a, a) , b, W are respectively as- 
sociated to the pi-electrons of each atom of the sublattices A and B and where the three 
relative vectors Si, <5 2 , S3 define the first nearest neighbors as depicted in figQ. These 
2D vectors are globally defined on the honeycomb and obey the remarkable constraint 
equation 5i + (5 2 -|-<53 = which, a priori, encodes also information on the electronic prop- 
erties of graphene. Throughout this study, we show amongst others, that the three above 
mentioned SU (2)'s are intimately related with these 5 n 's which, as we will see, are noth- 
ing but roots of SU (3). We also show that the mapping of the condition Yln=i $n = 
to the momentum space can be interpreted as a condition on the conservation of total 
momenta at each site of honeycomb whose solutions are classified by highest weight state 
representations of the SU (3) symmetry. We show moreover that the hamiltonian H of 
the tight binding model for first nearest neighbors has an interpretation in terms of the 
F ±/3 step operators of these SU (2)'s opening a window for more insight into the study 
of the electronic correlations in graphene and cousin systems. 

The organization of this paper is as follows: In section 2, we exhibit the SU (3) symme- 
try of graphene. In section 3, we give a field theoretic interpretation of the geometric 
constraint equation Si + S2 + £3 = both in real and reciprocal honeycomb. We also use 
the roots and weights of hidden SU (3) symmetry to study aspects of the electronic prop- 
erties of graphene. In section 4, we develop the relation between the energy dispersion 
and the hidden symmetries. In section 5, we give the conclusion and a perspective. 

2 Hidden symmetries of graphene 

In dealing with ideal 2D graphene, one notices the existence of a hidden SU (3) group 
symmetry underlying the crystallographic structure of the honeycomb lattice and govern- 
ing the hopping of the pi-electrons between the closed neighboring carbons. To exhibit 
this hidden SU (3) symmetry, let us start by examining some remarkable features on the 
graphene lattice and show how they are closely related to SU(3). 

Refereing to the two sublattices of the graphene monolayer by the usual letters A and B 
generated by the vectors a 12 = f(±v3, 3) and the relative ones 5i j2 = § (±V3, 1) , and 
denoting by 4> A (r^) and <p B (r,-) the wave functions of the corresponding pi-electrons, one 
notes that the interactions between the first nearest atoms involve two kinds of triva- 
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lent vertices capturing data on SU (3) symmetry, see fig© for illustration. This hidden 




Figure 2: (a) Nearest neighbors of a A- type atom, (b) Nearest neighbors of a B-type atom, 
(c) The fusion of the two vertices gives the adjoint representation of SU(3). 

SU (3) invariance can be made more explicit by remarking that the relative vectors d\, 
62 and 63 = —Si — 62 describing the three first closed neighbors to a A- type carbon 
at site Tj of the honeycomb together with their opposites — d n for B-type carbons are 
precisely the roots of the SU (3) algebra. Indeed, if forgetting about the scale dimension 
and thinking about the carbon-carbon distance a ~ 1.42 A° as the irrational numbeiQ 
a = a/2, one gets the remarkable identification 

61 = c*i , 8 2 = ol 2 j ^3 = -a 3 , (2.1) 

with Qi and ct2 being the two simple roots of SU (3). This symmetry can be also 
exhibited by computing the intersection matrix <5j • 8j of the two generators <5i and 62, 

(2.2) 

which turns out to be proportional to the Cartan matrix of the SU (3) algebra, Aij = 
a.i • ctj. Recall that the Lie algebra of SU(3) has rank two, eight generators {F a } and 
commutation relations that read in the Gell-Mann basis like 

[F a , F b ] = if abc F c , (2.3) 

with antisymmetric structure constants as /123 = 1, f 147 = f 516 = f 2m = f 257 = 
J345 _ y637 _ i anc j y458 _ j678 _ j__ j n ^ e Cartan- Weyl basis, useful for physical 

interpretations, we take h 1 = \/2F 8 and h 2 = \/2F 3 as the two Cartan terms and the 
six step operators like y/ZUf = F 1 ± iF 2 , ^/2U 2 ± = F 4 ± iF 5 , yfwf = F 6 ± iF 7 . The 
new commutation relations following from ( 12. 3 p read as follows, 

[h\Ut] = ±5\U± , (2.4) 




^Although it isn't a necessary condition for our analysis, this number is mysterious as it corresponds 
precisely to the length ||aj|| 2 = 2 of the roots of simply laced Lie algebras; see also end of conclusion. 
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with i — 1,2 and where the 5 n 's are the same vectors as in the graphene. Notice in 
passing that the hidden SU(3) described above seems to be just a sub-symmetry of a 
larger one since the <5 n 's obey the constraint relation 

3 

J2 S n = 0. (2.5) 



n=l 



which might hide an affine SU (3) Kac-Moody symmetry [HJ US] since the intersection 
matrix of the three relative 6i- vectors reads as <5j ■ Sj 



\Aij with 



/ 2 



A, 



\ 



-1 \ 
-1 

1 ) 



(2.6) 



describing exactly the generalized Cartan matrix of affine SU (3). Below, we shall restrict 
our study to the hidden ordinary SU (3) symmetry of the graphene and think about (12.51) 
as a physical constraint equation governing the electronic properties of the graphene. 



3 Electronic properties and SU (3) symmetry 

Quantum mechanically, there are two approaches to deal with the geometrical constraint 
relation ( 12. 5p . The first one is to work in real space and think about it as the conservation 
law of total space-time probability current densities at each site rj of the honeycomb. 
The second approach relies on moving to the dual space where this constraint relation 
and the induced electronic properties get a remarkable interpretation in terms of SU (3) 
representations. 



3.1 Conservation of total current density 

In the real space, the way we interpret eq (l2.5l) is in terms of the relation between the 
time variation of the probability density p(t, r*) = |0(£,rj)| 2 of the electron at site 
and the sum 

3 

^ t J Sn {t,T i )=J{t,T i ) , (3.1) 

71=1 

of incoming and outgoing probability current densities along the 8 n - directions. On one 
hand, because of the equiprobability in hopping from the carbon at to each one of the 
three nearest carbons at + 6 n , the norm of the Js n - vector current densities should be 
equal and so they should have the form 

J<5n (t, = j (t, rj) e n , n = 1,2,3 . (3.2) 
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These probability current densities together with the unit vectors e n = ^ point in the 
different S n - direction; but have the same non zero norm: HJ5JI = ||J,5 2 || = ||J<s 3 || = \ j\- 
Substituting in the above relation, the total probability current density J (t, r^) at the 
site r and time t takes then the factorized form 

JM = ^fe* B V (3-3) 

On the other hand, by using the Schrodinger equation ih^ = ^— ^ V 2 + V j <p describing 
the interacting dynamics of the electronic wave at r, we have the usual conservation 
equation, 

dp (t, r) 



dt 



+ divJ(t,r) = , (3.4) 



with probability density p (t, r) as before and J = ^p- (0V0* — 0* V0) with m the mass 
of the electron and <f> = <f> (t, r) its wave. Moreover, assuming §f = corresponding to 
stationary electronic waves <fi (t, r) = e luJt (j) (r), it follows that the space divergence of the 
total current density vanishes identically; div J = 0. This constraint equation shows that 
generally J should be a curl vector; but physical consideration indicates that we must 
have 

J(i,r) = , (3.5) 

in agreement with Gauss-Stokes theorem J y divJ dV = J dv J.d(T leading to the same 
conclusion. Combining the property J (t, r) = with its factorized expression | (J2 n 
given by eq f)3.3l) together with j 7^ 0, we end with the constraint relation J2 n $n = 0. 

3.2 Conservation of total phase 

In the dual space of the electronic wave of graphene, the constraint relation (12. 5ft may 
be interpreted in two different, but equivalent, ways; first in terms of the conservation of 
the total relative phase A(p tot = ^ k. Ar of the electronic waves induced by the hopping 
to the nearest neighbors. The second way is in terms of the conservation of the total 
momenta at each site of the honeycomb. 

Decomposing the wave function (r), associated with a A- type carbon at site r, in Fourier 
modes as ^ k e l2nV ' r 0(k); and similarly for the B-type neighboring ones <fi (r + 8 n ) = 
J2k e i27rk r (j) n (k) with k = (k x , k y ), we see that (f> (k) and the three 4> n (k) are related as 

(fc) = e i27r6 ^(k) , n = 1,2,3 , (3.6) 



with relative phases 9 n = k-(5 n . These electronic waves have the same module, 



but in general non zero phases; 6\ ^ 9 2 7^ #3. This means that in the hop of an 
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electron with momentum p = ftk from a site to the nearest at + 8 n , the electronic 
wave acquires an extra phase of an amount 8 n ; but the probability density at each site 
is invariant. Demanding the total relative phase to obey the natural condition, 

0! + 2 + 03 = , mod(27r), (3.7) 

one ends with the constraint eqf l2.5l) . Let us study two remarkable consequences of this 
special conservation law on the 9 n phases by help of the hidden SU (3) symmetry of 
graphene. 

(1) Using eq f)2.ip . which identifies the relatives 8 n vectors with the roots cx n of SU (3) 
symmetry, as well as the Lie algebra duality relation 

ct i -X j = 5 ij , i,j = 1,2 , (3.8) 

mapping the two simple roots oli, ct2 into the SU (3) fundamental weights Ai, A2, we can 
invert the three equations 9 n = k<5 n to get the p n =Kk n momenta of the electronic waves 
along the <5 n -directions. For the two first n 's, that is n — 1, 2, the inverted relations are 




Figure 3: Root /weight duality linking roots and weights CKj • dj = S^j. 



nicely obtained by decomposing the 2D wave vector k along the Ai and A2 directions of 
the dual lattice to end with the following particular solution, 

kx = 6 1 \ 1 , k 2 = 9 2 A 2 . (3.9) 

More general solutions of type ki = 6\ X\ + x 2 A2 and k 2 = X\ X\ + 9 2 A2 will be consid- 
ered in next subsection. Notice by the way that the 2D vectors Ai and A2, interpreted 
in the framework of SU(3) group theory as the fundamental weights, are nothing but 

A 1 = | ai , A 2 = |a 2 , (3.10) 

where a 12 = |(±v3, 3) stand for the generators of the A- type atoms introduced in the 
beginning of section 2; see also fig©. The above relations 1 13. 10p give an other evidence 
for the role of the SU(3) symmetry in the study of the electronic properties of graphene. 
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(2) To get the wave vector from the relation 83 = k ■ 63, we decompose the 2D vector 
like k = q 1 X 1 + q 2 A 2 ; then substitute <5 3 = —oci — et 2 and use eq (l3.8p to end with 
#3 = —qi — q%. Comparing with eq (!3.7[) . we find that the wave vector k 3 of the electronic 
wave along the £3- direction reads as follows, 

k 3 = 6>xAi + 6> 2 A 2 . (3.11) 

Now, combining (13.91) and (13.111) . we find that the property 9„ = describing the 
conservation law (13. 7\\ of the total phase of the electron hops to nearest neighbors can 
be mapped to a constraint relation on the conservation of total outgoing and incoming 
momenta Kk n at each site of the honeycomb, i.e: 

k 1 + k 2 -k 3 = . (3.12) 

This result is not strange; it may be directly obtained by mapping (12.51) to the reciprocal 
lattice. Below, we study the solutions of this constraint relation in connection with the 
hidden SU (3) symmetry of the honeycomb. 



3.3 More on the constraint eq (!3.12j ) 

Seen that the relative vectors 8 n defining the first nearest neighbors are roots of SU (3) 
as shown by eq (l2.1l) , a way to deal with the constraint relation (I3.12j) is to think about it 
as a SU (3) group representation relation. This means that the wave vector in eq (l3.11l) 
may be thought of as given by kx = £ x Ai, k 2 = £ 2 A 2 , k 3 = £ 3 A 3 where the £ n 's are 
real numbers and A 1; A 2 , A 3 are three generic weight vectors of SU(3). An interesting 
situation corresponds to the case where ^ = £ 2 = £ 3 = £ allowing to turn the constraint 
eq (l3.12l) into a constraint relation on SU (3) weights, 

£(A 1 + A 2 -A 3 ) = . (3.13) 

This is a remarkable relation which may be motivated by thinking about SU (3) as a 
basic symmetry that governs the electronic properties in graphene. After all, eq( 13.13l) is 
the dual of (12.51) and moreover Xm=i <5n = is itself a SU (3) condition; see also footnote 
1. Under this hypothesis, and thinking about the A n 's as highest weight vectors that 
can be decomposed as, 

A„ =p„iAi +p„ 2 A 2 , (3.14) 

we get, after substituting in (13.131) . the following conditions on the pij positive Dynkin 
integers, 

£(pn +P21 -P31) = , ^ 15 « 

£ (P12 +P22 -P32) = . 
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The simplest solution of these relations corresponds to taking the weight vectors as 



Ax = Ai , A 2 = A 2 , A 3 = X adj = Ai + A 2 , (3.16) 

in agreement with eqs fl3.9H3.lip . Particular solutions type A = pXi and A' = p\ 2 with 
positive integer p are in the same class as Ai and A 2 . From this analysis we learn that 
within the SU (3) set up, the solutions of ( 13. 12ft have the following features: 

(i) the norm of the wave vector of the (p (k n ) wave is = | (p^ + p 2 n2 + p n \p n 2) £ 2 , 

(ii) the phases of the waves 4> n (k x , k y ), defined by tp n = arctan (j^j, are quantized as 

knx Pnl~Pn2 \ ' J 

(iii) the Fourier waves (f> (k x , k y ) may be interpreted as ID field multiplets transforming 
into SU (3) highest weight representations as, 

0(^A) = 0a(O = (m) (O , (3- 18 ) 

with dimension 

l(l+p)(l + g )(2+p + g) . (3.19) 

In this SU (3) picture, the physics describing the electron hops between the nearest 
carbons is completely captured by SU (3) highest state representations. For instance, 
taking Ai = 3 and A 2 = 3 as in figfjlj), and using the tensor product decomposition 

3®3 = 9 = 8©1 , (3.20) 

the Fourier waves propagating between the nearest sites are 0(i o)> 0(0,1) > 0(i,i) anc ^ ^(o,o)- 
The zero mode ^mo) describes the state where the electron doesn't hop. Notice that 
each of these A (£) fields hide sub- modes (4>\) fl associated with the various states of the 
representations. For the example of fig®, we have: 



^A = 0A 



iX—ai 


<t>\ 










,±(Vi-ai, 2 ) 






0i° 



(3.21) 



with X adj = ip. This de generacy can be lifted by breaking down the hidden SU (3) sym- 
metry of the graphene. This may be achieved by implementing defects in the honeycomb 
that violate eq (l3.15|) . 
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Figure 4: Real and dual lattices represented simultaneously. Total probability current density 
and total momenta of the electronic waves are conserved at each site of the honeycomb. 



4 Energy dispersion and the hidden symmetries 

The energy hi describing the hopping of a pi-electron from the site to its three nearest 
neighbors at rj + S n is nicely represented by the tight binding hamiltonian [16] whose 
total form reads as H = — t ^ hi, 

2 

hi = ^2 o. ( r fot ( r i + On) + he, (4.1) 

n=0 

where cij, bj, a\, bj are fermionic annihilation and creation oscillators and t ~ 2.8eV the 
hopping energy. With this hamiltonian H, one learns much about the electronic band 
structure of graphene. However to get more insight about the hidden symmetries of the 
honeycomb, it is interesting to express H in terms of the i^ ±a ™ steps operators gener- 
ating SU(2) sub-symmetries inside SU(3). To do so, we start from the wave functions 
cj) A (r^ =< r\(j) i > and 4> B (r« + a n ) =< r|0f n > associated with a fixed A-type atom 
and its nearest B- type neighbors. Then use the structure of the honeycomb (figdTJ)) to 
write down the action of the F i a "'s generating the electron hopping. At each of the 
sublattice A, we have 



K an \^) =\€ n ) , i^ +Q "IC") = o , 

Fr an \€ n ) = l&> , Fr- lb) = , 
from which we read the following relations, 

[^iT 1 ] =K n , [hr,Ft a -] = ±Ft a " , 

{Ft an ,Ff an } =0 , {F<* n ,Fr a "} = J«™ . K ' } 
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with J" n a commuting central element [i^ " n , J" n ] = 0. The commutation relations tell 
us that locally each set (/if", F i an ) generate an SU (2) group along the a n -direction 
in the hidden SU(3) symmetry. The anti-commutation relations, which read also like 
(F^"™) 2 = 0, requires F[° an to be in the isospin | representation; that is 2 x 2 matrices 
linking the two sublattices A and B of the honeycomb. The fermionic realization (14.11) 
is a representation of eqs fl4.3l) where the i^ ±a "'s are solved as 

Ft = (lib} , hi = bib} - did] 

F i = a % > Ji = b M + a i a l 

In terms of the globally defined operators F ±OLn = '^2 i Ft° Ln , the hamiltonian H takes 
the simple form, 

2 

H = -t i F+a " + F ~ an ) > ( 4 - 5 ) 

n=0 

where we have used a = —oc 3 . Besides hermiticity, H has two special features that 
we want comment: (i) H is not invariant under the SU (2) symmetries since along with 
(14. 5 p we also have the cousin operators 

2 2 

L = -itfJ2 ( F+an ~ F ~ an ), M = E h an , ( 4 - 6 ) 

n=0 n=0 

obeying the commutation relations 

[fJl^M, [M,f]=& [M,§]=f. (4.7) 

The real number $ in (14. 6p may be interpreted in terms of coupling to a constant external 
magnetic field, (ii) H is not a positive definite operator in the sense that its energy 
spectrum has two signs; a region with positive energy describing the conduction band 
and a region with negative energy associated with holes. 

Performing the Fourier transform of the step operators F +OLn = ^ e l2 n G k and putting 
back into the hamiltonian, we can put H in various forms; in particular like 

H = -tJ2(i> k Gk + ^ k G k ), (4.8) 

with ifj k = e i2nkai + e i27rk '" 2 + e -*&*-«3. Setting Q = e i27rf and k = £A with A = 
PiAi + P2X2, we can bring this hamiltonian to —t A H\ (£) with 

H A = ijf 3) G A + ijf 3) G A , (4.9) 

and i{j a {3) = [Q A - ai + Q Aa2 + Q-A.oi-A.a a j _ To get the wave vectors k f = £ F A F at 
the Fermi level, one has to solve the zero energy condition ifj^ 3 ^ (£) = whose solutions 
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(4.4) 



are given by the cubic root of unity (1 + Q + Q 2 ) = 0. They are generated by £ F = | 
times the fundamental weights of the dual lattice; i.e = |Ai and = modulo 
translations. Moreover setting 

G k = e* r °V k , V fc = l^fc|e w * , (4.10) 
where Vk is positive definite and Tq hermitian; then substituting in ( 14.5}) . we get 

h = J2 £ kN k - £kMk > ( 4 - n ) 



with Ek = \t^k\ anc ^ 



M fc = 2cos 2 (^AV, 



k = z cus- \ — -= — 

- V (4.12) 

fc 



A4 = 2sin 2 (**±^)p, 

We end this section by first noting that using the fermionic realization, eq (14.12j) gets 
a simple interpretation in terms of electron and hole number operators ej.e/. and h k h k - 
Second, the group theoretical approach developed in this study may be used to deal with 
graphene multilayers. In the case of graphene bilayer, one expects symmetries of type 
SU (3) x SU (2) x SU (3) with each SU (3) factor as before and where the SU (2) term 
refers to transitions between the two layers. 



5 Conclusion 

In this paper, we have shown that 2D graphene has a remarkable hidden SU (3) symmetry 
that allow to classify the propagating Fourier waves (f>(k x ,k y ) in terms of ID highest 
weight field multiplets (p A (£). Conservation of total incoming and outgoing momenta 
at each lattice site translates into triplets of SU (3) HWRs constrained by eq (l3.15[) and 
whose basic one is (3,3,8) satisfying 3 <g> 3 = 8 © 1. We have also shown that, from 
SU (3) view, tamp = jN- is quantized as a/3 and moreover the hamiltonian of the 
tight binding model is the sum of the step operators of the three SU (2) sub-symmetries of 
the hidden SU (3) invariance. This connection with Lie algebra teaches us that graphene 
may be thought of as the second element of class of theoretical systems involving higher 
rank symmetries [151 EE]- The first element has a hidden SU (2) and energy spectrum 
£ a (0 = i^/^V'A 4 i this should correspond to the linear poly-acetylene chain with 
hamiltonian H su ^ — — t hf 1 ^ and, 

hT {2) = a (r<) 6 f (r< + a) + a (r<) 6 f (r< - a) + /ic, (5.1) 

where a stands for the S 1 ?/ (2) root. It is also interesting to note the connection between 
Yln=i $n — and the imaginary root of affine SU (3) suggesting that ideal graphene 
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could have a richer hidden symmetry containing SU (3) as the zero mode. If this is 
the case, graphene would also exhibit a hidden 2D conformal structure underlying the 
honeycomb lattice and capturing information on eventual critical behaviors of electronic 
correlations along the line of (TTJ [15]. This issue is understudy; it will be developed 
elsewhere. 
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